Simulating Cardiac Fluid Dynamics in the Human Heart.

Cardiac fluid dynamics fundamentally involves interactions between complex blood flows and the structural deformations of the muscular heart walls and the thin, flexible valve leaflets. There has been longstanding scientific, engineering, and medical interest in creating mathematical models of the heart that capture, explain, and predict these fluid-structure interactions. However, existing computational models that account for interactions among the blood, the actively contracting myocardium, and the cardiac valves are limited in their abilities to predict valve performance, resolve fine-scale flow features, or use realistic descriptions of tissue biomechanics. Here we introduce and benchmark a comprehensive mathematical model of cardiac fluid dynamics in the human heart. A unique feature of our model is that it incorporates biomechanically detailed descriptions of all major cardiac structures that are calibrated using tensile tests of human tissue specimens to reflect the heart's microstructure. Further, it is the first fluid-structure interaction model of the heart that provides anatomically and physiologically detailed representations of all four cardiac valves. We demonstrate that this integrative model generates physiologic dynamics, including realistic pressure-volume loops that automatically capture isovolumetric contraction and relaxation, and predicts fine-scale flow features. None of these outputs are prescribed; instead, they emerge from interactions within our comprehensive description of cardiac physiology. Such models can serve as tools for predicting the impacts of medical devices or clinical interventions. They also can serve as platforms for mechanistic studies of cardiac pathophysiology and dysfunction, including congenital defects, cardiomyopathies, and heart failure, that are difficult or impossible to perform in patients.

A model of fluid-structure and biochemical interactions for applications to subclinical leaflet thrombosis.

Subclinical leaflet thrombosis (SLT) is a potentially serious complication of aortic valve replacement with a bioprosthetic valve in which blood clots form on the replacement valve. SLT is associated with increased risk of transient ischemic attacks and strokes and can progress to clinical leaflet thrombosis. SLT following aortic valve replacement also may be related to subsequent structural valve deterioration, which can impair the durability of the valve replacement. Because of the difficulty in clinical imaging of SLT, models are needed to determine the mechanisms of SLT and could eventually predict which patients will develop SLT. To this end, we develop methods to simulate leaflet thrombosis that combine fluid-structure interaction and a simplified thrombosis model that allows for deposition along the moving leaflets. Additionally, this model can be adapted to model deposition or absorption along other moving boundaries. We present convergence results and quantify the model's ability to realize changes in valve opening and pressures. These new approaches are an important advancement in our tools for modeling thrombosis because they incorporate both adhesion to the surface of the moving leaflets and feedback to the fluid-structure interaction.

Patient-Specific Immersed Finite Element-Difference Model of Transcatheter Aortic Valve Replacement.

Transcatheter aortic valve replacement (TAVR) first received FDA approval for high-risk surgical patients in 2011 and has been approved for low-risk surgical patients since 2019. It is now the most common type of aortic valve replacement, and its use continues to accelerate. Computer modeling and simulation (CM&S) is a tool to aid in TAVR device design, regulatory approval, and indication in patient-specific care. This study introduces a computational fluid-structure interaction (FSI) model of TAVR with Medtronic's CoreValve Evolut R device using the immersed finite element-difference (IFED) method. We perform dynamic simulations of crimping and deployment of the Evolut R, as well as device behavior across the cardiac cycle in a patient-specific aortic root anatomy reconstructed from computed tomography (CT) image data. These IFED simulations, which incorporate biomechanics models fit to experimental tensile test data, automatically capture the contact within the device and between the self-expanding stent and native anatomy. Further, we apply realistic driving and loading conditions based on clinical measurements of human ventricular and aortic pressures and flow rates to demonstrate that our Evolut R model supports a physiological diastolic pressure load and provides informative clinical performance predictions.

Failure properties and microstructure of healthy and aneurysmatic human thoracic aortas subjected to uniaxial extension with a focus on the media.

Current clinical practice for aneurysmatic interventions is often based on the maximum diameter of the vessel and/or on the growth rate, although rupture can occur at any diameter and growth rate, leading to fatality. For 27 medial samples obtained from 12 non-aneurysmatic (control) and 9 aneurysmatic human descending thoracic aortas we examined: the mechanical responses up to rupture using uniaxial extension tests of circumferential and longitudinal specimens; the structure of these tissues using second-harmonic imaging and histology, in particular, the content proportions of collagen, elastic fibers and smooth muscle cells in the media. It was found that the mean failure stresses were higher in the circumferential directions (Control-C 1474kPa; Aneurysmatic-C 1446kPa), than in the longitudinal directions (Aneurysmatic-L 735kPa; Control-L 579kPa). This trend was the opposite to that observed for the mean collagen fiber directions measured from the loading axis (Control-L > Aneurysmatic-L > Aneurysmatic-C > Control-C), thus suggesting that the trend in the failure stress can in part be attributed to the collagen architecture. The difference in the mean values of the out-of-plane dispersion in the radial/longitudinal plane between the control and aneurysmatic groups was significant. The difference in the mean values of the mean fiber angle from the circumferential direction was also significantly different between the two groups. Most specimens showed delamination zones near the ruptured region in addition to ruptured collagen and elastic fibers. This study provides a basis for further studies on the microstructure and the uniaxial failure properties of (aneurysmatic) arterial walls towards realistic modeling and prediction of tissue failure. STATEMENT OF SIGNIFICANCE: A data set relating uniaxial failure properties to the microstructure of non-aneurysmatic and aneurysmatic human thoracic aortic medias under uniaxial extension tests is presented for the first time. It was found that the mean failure stresses were higher in the circumferential directions, than in the longitudinal directions. The general trend for the failure stresses was Control-C > Aneurysmatic-C > Aneurysmatic-L > Control-L, which was the opposite of that observed for the mean collagen fiber direction relative to the loading axis (Control-L > Aneurysmatic-L > Aneurysmatic-C > Control-C) suggesting that the trend in the failure stress can in part be attributed to the collagen architecture. This study provides a first step towards more realistic modeling and prediction of tissue failure.